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Collective charge transport is studied in one- and two-dimensional arrays of small 
normal-metal dots separated by tunnel barriers. At temperatures well below the 
charging energy of a dot, disorder leads to a threshold for conduction which grows 
linearly with the size of the array. For short-ranged interactions, one of the correlation 
length exponents near threshold is found from a novel argument based on interface 
growth. The dynamical exponent for the current above threshold is also predicted 
analytically, and the requirements for its experimental observation are described. 
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Systems exemplifying collective transport in quenched disorder include sliding charge- 
density waves (CDW's) fluids in disordered media ||, and type-II superconductors 
For these dynamical systems, there does not yet exist a classification whereby the long- 
wavelength behavior can be predicted from the characteristics of the microscopic degrees of 
freedom. To study this question of universality experimentally requires systems where the 
microscopic degrees of freedom, the range of interactions, and the nature of the disorder 
are well understood. Here, we propose as a model system an array of small-capacitance 
normal-metal dots. 

In this letter we examine the low-temperature, nonlinear charge transport in such an 
array. The dots are treated as capacitively coupled conductors with charges allowed to 
tunnel between neighboring dots. In contrast with previous work |J, we explicitly include 
the effects of random offset charges on each dot and investigate the limit where the number of 
dots becomes large. We find that the onset of conduction occurs at a voltage Vr proportional 
to the linear array size. One of the correlation lengths that diverges near this threshold 
is found from a general argument based on interface growth, while another is found by 
focusing on "slow points" which control the current. These correlation lengths determine 
the branching of current paths in the array and hence the current near onset. In particular, 
we predict that the current through linear and square arrays behaves as 

/ ~ (y/V T - l) c (1) 

with ( = 1, 5/3 in dimensions d — 1, 2, respectively. 

The array we study is depicted in Fig. |l|. For a tunneling resistance R between dots 
large compared to the quantum resistance h/e 2 , the state of the array is fully described by 
the number of electrons in each dot. The energy is then all electrostatic and is determined 
by a matrix of capacitances CV,-. We assume a constant capacitance C between neighboring 
dots and between the leads and adjacent dots, and a capacitance C g between each dot 
and the back-gate which underlies the entire array. The leads and back-gate are taken to 
have infinite self-capacitance. We concentrate on the Coulomb-blockade regime, where the 
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thermal energy is much smaller than the charging energies, i.e., ksT <C e 2 / [2max(C, C g )} 
(for 1/im dots, this energy is ~ lmeV ||). We measure distances in units of the dot spacing. 
Given the charge Qi on each dot, the electrostatic energy is 

E = E i.Qi + <l>) ( ',/( ( li + Qj) + VlQl + V R Q R + £ (Va + V iR - V g )Q u (2) 

dots i,j dots i 

where Ql,r are the charges on the leads, at voltages V^r, and V g is the back-gate voltage. 
Disorder is included through the offset charges which represent the effective charge on each 
dot due to nearby charged impurities. Large fluctuations in disorder will be compensated 
by an integral number of mobile charges, so that < qi < e. The voltage on dot i due to 
the left (right) lead is given by Vh(r) = Cojr\Vl(R)- In general, the elements of the inverse 
capacitance matrix fall off exponentially with a screening length A that increases with C/C g 
(for C » C g , A « (C/C,) 1 / 2 ). 

At low temperatures, a charge may tunnel between dots only if such an event lowers the 
electrostatic energy of the array. The kinetic energy gained by the tunneling electron is as- 
sumed to be dissipated [0. The tunneling rate from one configuration S = (..., Qi, Qj, ...) 
to another configuration S' = (■■■,Qi — 1, Qj + 1, ...), where i and j are neighboring dots 
or a dot and a neighboring lead is given by 

v s ^s> = (e'R)- 1 9 (E(S) - E(S')) [E(S) - E(S')} . (3) 

This rate grows linearly with energy gain since the number of electrons available to tunnel 
is proportional to the relative shift of the Fermi surfaces in the dots. For arrays of a few 
junctions, numerical results on this model compare well with experiment 0. 

For large arrays, we find a second-order transition, with associated critical phenomena, 
which separates a static, non-conducting state from a dynamic, conducting state. The 
control parameter is the voltage difference between the leads. At low voltage differences, 
the array always relaxes to a static configuration, while at high voltage differences, charges 
traverse the array from one lead to the other. 

An important question to ask is whether the conduction transition is hysteretic for a 
given realization of disorder. In one-dimensional systems, the current is a unique function 
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of the applied voltages, regardless of the magnitude of A ||. In two-dimensional arrays at 
zero temperature, the current can depend on the history of the applied voltages. To within 
our numerical accuracy, however, the current is history independent in typical samples. It 
can furthermore be shown that the current is entirely independent of history in the limit of 
short screening lengths, C/C g — > 0, to which we shall devote most of our attention. 

We have numerically determined the dependence of the threshold voltage for conduction, 
Vt(N) = Vl — Vr, on the ratio C/C g and on the linear system size N (for fixed gate voltage 
V g = 0). We find that the threshold voltage is proportional to N, 

lim MN)C g /Ne = a(C/C g ), (4) 

iV — >oo 

where the overbar represents an average over disorder. 

The function a(C/C g ) for one-dimensional arrays is plotted in Fig 0. In the limit C/C g — > 
0, the voltage on a dot is just (Qi + qi)/C g , as the capacitive coupling between dots is 
negligible. The schematic in Fig. 0(a) shows that in order to carry a current in this limit, 
the voltage difference across the array must be large enough to overcome ~ N/2 upward 
steps in the random potential. This observation gives a(C/C g — >• 0) = 1/2. In the limit 
of large C/C g (large A), Vp can be estimated by balancing the "force" on the charges due 
to a charge density gradient against the random potential gradient f|. It is necessary 
to recognize that there is a stability limit for the dot-to-dot potential difference: at higher 
potential differences, charges will tunnel, reducing the potential across the tunneling barrier. 
Estimating the magnitude of the pinning forces to be given by this stability limit, we find 
T0|| a maximum static density gradient of ~ e/A 2 , that is, a density change of 0(1) charge 



per screening length is allowed in regions separated by A. This gives a ~ (C/C g )~ l at large 
A. Numerically, we find in d = 1 that (C/C g )a — ► 0.10(1), as C/C g — > oo. 

We now discuss the approach to the conduction threshold in two dimensions, in order 
to elucidate the critical behavior of the correlation lengths and current. We concentrate on 
the limit C/C g — > 0, that is, A small, in order to maximize the number of effective degrees 
of freedom and to eliminate hysteresis. Taking Vr to be fixed and raising Vl, charge moves 



from the left lead onto the array. The condition for a charge to overcome the Coulomb 
barrier and tunnel from site i to neighboring site j is 

V t > V, + e/C 9 . (5) 

At a given Vl, this advance of charge is halted when Vi < Vj + e/C g everywhere. Though 
the tunneling is stochastic, the static configuration at any Vl < Vt is entirely determined by 
the disorder realization. The distance to which charge penetrates therefore defines a unique 
interface (given Qi = initially). When Vl is raised by e/C g , charge is added to each point 
on the interface and therefore the interface must advance by at least one lattice spacing. In 
addition, the interface may advance further at some points if the local disorder is favorable. 
The motion of this interface is depicted in Fig. 0(b); conduction occurs when the interface 
reaches the right lead. Numerical calculations of this threshold give a(C/C g — 0) = 0.338(1) 
in al = 2. 

The dynamics of the Coulomb-blockade condition Eq. (||) make the interface "motion" 
with increasing Vl similar to the stochastic growth of interfaces in models without quenched 
spatial disorder, such as the one due to Eden [IT|. The results on the Kardar-Parisi- Zhang 



(KPZ) equation for a (d — l)-dimensional interface |T2j subject to short-range correlated 
noise are therefore useful in understanding the behavior of a <i-dimensional array of dots. 
This is to be contrasted with the usual motion of interfaces at small velocities through 
random media, where the interface can be pinned for some time at one point, resulting in 



long time correlations and exponents distinct from those of the KPZ equation |T3 . 

In the case d = 2, the results for KPZ interfaces imply that the width of the interface 
must scale as V L ||12|| . Furthermore, the fluctuations in the position of maximum advance 



of the interface behave as ~ Vl r \\nV L ) l/2 . The rms fluctuations AVr in Vt and the mean 
threshold voltages Vt therefore behave as 

AVt/W ~ iV- 2 / 3 (ln iV) 1/2 , (6) 
Vr(N)C g /eN -a~ N- 2/3 (\nN) 1/2 } (7) 
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for an N x iV array flj (for d = 1, AVt/Vt ~ A^' 1 / 2 ). The fluctuations in the threshold 
voltage as a function of size may be used to define a finite-size scaling exponent via 
N ~ (AVt/Vt) - ^- This length, besides giving the fluctuations in Vr, determines the finite 
size crossover in quantities such as the polarization of the array ||. From Eq. @, we find 
vt — 3/2. Numerical simulations on systems up to size iV = 2560 are fit very well by Eq. 
(|), as shown in the inset in Fig. 0(a), and by Eq. (|7|). 

The current in the one-dimensional model with only on-site interactions can be under- 
stood in detail. For voltages much greater than threshold, v = (Vl — Vt) /Vr 1, the 
charge gradient across each junction is much greater than one. By Eq. (|3p, the current is 
then approximately 

/ « (e/2RC g )v. (8) 

In contrast, near threshold, the discreteness of the charges and the disorder become impor- 
tant. The excess charge gradient above the threshold configuration is composed of steps 
that occur at well separated "slow points", located where the potential drop between dots 
in the threshold configuration is small compared to e/C g . It can be shown [IU| that the cur- 



rent is given by the fastest "slow point" ; the tunneling rate across this point is, on average, 
(V — Vr)/eRN. Interestingly, near threshold this also gives Eq. (g|). We therefore find ( = 1 
for Eq. (P in d — 1. As shown in Fig. 0(a), Eq. (§) is consistent with our numerical results 
near and far from threshold. 

The pattern of current flow in a typical two-dimensional array is shown in Fig. 0(b). At 
voltages just above threshold, Vl — Vt « e/C g , the current is in general carried on a single 
path, with little or no branching. This path is exactly the one with the minimal number of 
upward steps in the potential between the two leads. Previous work [jl2| shows that such 
paths have transverse fluctuations ~ n 2 ^ 3 , where n is the distance from the left lead; this 
is consistent with our numerical results. Increasing the voltage to a few times e/C g above 
threshold opens multiple channels which branch and reconnect, as shown in Fig. |3](b). Note 
that at voltages Vl exceeding threshold by 0(eA^ 1//3 (ln N) 1 ^ 2 /C g ), current can in principle 
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flow anywhere in the array, as the charge invasion interface contacts the right lead at each 
point. However, near threshold, all but a small fraction of the current is confined to a few 
major current-carrying paths. The selection of these paths out of all possible paths results 
from a characteristic length for branching, as we now discuss. 

Any current- carrying channel between the two leads must have (Vj, — Vr)C g excess steps 
in the charge density relative to the threshold configuration. This gives a correlation length 
£|| = eN/iVi — Vr)C g which separates the steps in the excess charge, as in d = 1. This 
length determines the separation between branch points along the channels fH| . The chan- 



2/3 

nels therefore wander transversely a distance £m between branch points, giving a channel 
separation of £j_ ~ v~ 2 ^ 3 . The current through each channel behaves as ~ (e/2RC g )v, since 
each segment between branch points is one-dimensional. The current through the array is 
then given by 

I ~ (e/2RC g )vN/£ ± ~ (e/2RC g )Nv 5/3 J (9) 

resulting in ((d = 2) = 5/3. 

Our numerical results for the transport in two-dimensional systems are shown in Fig. 
0(b). The current- voltage relationship is approximately fit by I/N ~ t> 2 ' , over the range 
10~ 2 < v < 10 _1 , but the slope on a log-log plot does not converge in the range we have 
studied numerically. To observe the true exponent requires arrays larger than 400 2 , either 
numerical or experimental. 

In conclusion, we have determined the threshold for conduction in arrays of small normal- 
metal dots with disorder. By examining correlation lengths that describe the separation of 
parallel current paths and the distance between dynamically important "slow points" we 
have determined the transport behavior near the threshold. The critical exponents for the 
current and correlation lengths which we have derived using the KPZ interface model are 
distinct from those found for elastic media |1|] and fluid flow |3|,[nj. These differences are 



clearly related to the novel features of this system, namely the discreteness of the carriers 
and (quantum) stochastic flow, which result in (a) an always advancing charge interface 



below threshold and (b) the non-local selection of current paths above threshold. 

We wish to thank Paul McEuen for encouraging our interest in arrays, and Chao Tang 
for many valuable discussions. 
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FIGURES 

FIG. 1. The threshold voltage per dot, Vt/N, in units of e/C g , for conduction through 
a one-dimensional array of normal- metal dots as a function of C/C g . The dashed lines show 
analytical predictions. A two-dimensional array of dots is shown in the inset. The indicated 
capacitance between dots is C, while the capacitance between each dot and the back-gate (the 
dashed rectangle) is C g . At T = 0, charges may only tunnel between neighboring dots if this 
lowers the total electrostatic energy. The voltages applied to the left lead, the right lead, and the 
gate are indicated as Vl, Vr, and Vq, respectively. 

FIG. 2. (a) Schematic of dot voltages for a one-dimensional array below the threshold for 
conduction, in the limit of short screening length A. Each square indicates an increase in the 
on-site voltage by e/C g due to an added charge; the relative offset in the voltages is caused by 
quenched disorder. If the left-lead voltage Vl is further raised by e/C g , charges will tunnel onto the 
array until stopped by the next upward step, as indicated by the dashed squares. The inset shows 
the calculated sample-to-sample fluctuation of threshold voltages, Vr, in one and two dimensions 
as a function of linear system size, N, with fits described in the text ( Eq. @). (b) Contours 
of constant charge occupation in a 160 2 array at the threshold Vp (contour spacing is 5 charges). 
Successive contours coincide with the distance to which charge flows for various voltages below 
threshold, Vl <Vt- Conduction occurs at the voltage where the charge first reaches the right lead. 

FIG. 3. (a) Plot of current-voltage relationship near threshold for one- and two-dimensional 
arrays of various sizes. The numbers in parentheses give the number of disorder realizations. For 
one-dimensional arrays, the current both near and far from threshold is well fit by Eq. (|8|). The data 
for the d = 2 arrays are approximately fit by I ~ (V — Vr) 2 ' at the lowest currents shown, though 
the local slope on the log-log plot has not converged, (b) Current paths in a two-dimensional array 
(of size 160 2 ) at two voltages near threshold. Very near threshold, the current flows in a single 
narrow channel (dark line). Multiple, branching channels are shown for a voltage 3e/C g above 
threshold (light lines). 
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Fig. 1 ~ Middleton and Wingreen, "Collective transport in arrays of quantumdots" 




Fig. 2 ~ Middleton and Wingreen, "Collective transport in arrays of quantumdots" 
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Fig. 3 — Middleton and WIngreen, "Collective transport in arrays of quantumdots" 



